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Abstract 

We discuss the low-temperature specific heat of the integrable SU(N)- invariant 
Heisenberg model in one dimension with degrees of freedom in the symmetric rank- 
in tensor representation, especially for the antiferromagnetic coupling. It is known 
that the linear specific heat coefficient function of a field which breaks 

the SU(N) invariance of the internal degrees of freedom. We calculate the ■y[ N ' m ~i in 
a zero field and in a small field. The in-field jl N > m i is less than the zero-field ^ N ' m ^ 
as expected since the entropy is reduced in the ordered system. The zero-field 
ryW,m] j g same as tn e one obtained by the prediction of the critical behavior 
of ID quantum spin systems via conformal field theory. This extends the previous 
results for N=2 to an arbitrary N. 



The Bethe-ansatz method yields singular behavior in numerous other models 
with SU(N) symmetry, for instance, the N-component ID fermion gas interacting via a 5- 
functional repulsion the Gross- Neveu model 0, the N-component supersymmetric 

ID t-j model |J, and the N-degenerate ID Hubbard chain [0, §]. In particular, the 
logarithmic singularities in zero-temperature susceptibility to a field are common to these 
models. For a small Zeeman splitting H and for arbitrary N and m j9|, 
x(H) ( m m 2 ln\NlnH\ \ 

X(0) = V + N\lnH\ NVH + " ' J ' ^ 
The nonanalytic field dependence of the susceptibility suggests a singular behavior in the 



low-temperature specific heat, C = 7T, as H and T tend to zero P]-|TT|, i.e. 



lim lim 7 [7V ' m] ^ lim lim 7 [Ar ' m] . (2) 

In this letter we calculate the linear specific heat coefficient yt N > m i in the limit 
of T — > with and without the SU(N)-invariance breaking field. Our main result is that 
for a small magnetic field H, 



limlim7M = 

H^OT^O 12 J 



1 + 



N- 1 



, , , rv™i N(N 2 - l)m , s 

^ lim lim y [N ' m] = — '—, 3 

^ t^oh^o 1 6J(m + N) V ; 



where 



Mm, j)(l)) + o V^ J KN - j) r(i - ffi)r(i - 



2-k' \ 2n 

2nej m (JV — 

Q>N,m,j at ) 

with T being the gamma function. 

The generalized SU(N) Heisenberg model in one dimension with degrees of free- 



dom in the symmetric rank-m tensor representation is given by [12|, [T3 

N a 



H^ = JJ2l^l (5) 

k=l 



with 

m j—1 m 

&i = EE^-^r 1 II <(A fe ®A fe+1 + /f H ), 
0=1 s=0 i=o,i^j 

pl N ' m] = l(2m+l-l)-m 2 (N -1)/N, 

where J > is the antiferromagnetic coupling constant, k labels the N a lattice sites in 
the chain and A fc denotes the generators of the symmetric rank-m tensor representation, 
playing the role of spin operators. 

The Bethe's ansatz method applied to this model yields N — 1 sets of rapidities 
and all rapidities within a given set are different. Hence rapidities obey Fermi statistics. 
In general, these rapidities are complex and form so called, strings for a very large system 
such that the center of mass of each string becomes real. At finite T the length n of strings 
could be infinitely long. In the thermodynamic limit we define the densities of states for 
occupied rapidities as ay in (A) and for unoccupied ones as <r^ n (A). These densities then 
satisfy the following integral equations with n — 1, • • • , oo, 

N-l 

*r,n(A) + <(A) - E PrZ * «+l + ^n-l " ^-l,n " *s + l,n) = (X)5 nm , (6) 
s=l 

where ao,n = CN,n = for all n, the asterisk denotes a convolution, P^ s is the Fourier 
transform of 

~ N sinh(mm(r, s)^) sinh([iV — max(r, s)]f) 

r ' s ^ = sinh(f) sinh(TVf) ' ^ 



and 



1 sin( it - N ~ r ' 



P*r (^) at i/2ir\\ / N—r\ " (^) 

iv cosh(^A) + COS^^^) 

It is usual to introduce the thermodynamic energy potentials £ r , n (A) such that 

<J rn(^)/°"r,n(A) = exp{e r , n {\) /T) . Hence the strings with the same lengths in each family 



3 



form energy bands denoted by £ r , n (A). These thermodynamic energy potentials are de- 
termined by the following infintely and non- linearly coupled integral equations, so called, 
the thermodynamic Bethe ansatz equations, 

N ~ 1 fl + p £.,»+i/TWl I p e s ,n-i/T\ 

£r ,„(A) + T £ P« * In 1^,,..,^;^,^ = -l*JK»(\)6 n , m , (9) 

where £o,n = £N,n = — oo for all n. It is important to notice that Eq. (|6|) and Eq. (|9|) 
consist of the same integral kernels and driving terms. This property makes the zero-field 
specific heat be obtained analytically at sufficiently low temperature. We will discuss it 
in details later. 



At T=0 Johannesson |L3| have shown that only strings of length m contribute. 
Note that m is nothing but the rank of the tensor representation of the model. In other 
words we have only one band (n = m) for each family (r) in the ground state. Hence each 
band is described by the thermodynamic energy potential £ rim (A) which is symmetric and 
monotonically increasing function of A. e rjTn (X) changes sign when symmetry breaking 
field H r is applied. This means that in field £ r , m (A) crosses zero at a special value of A, 
designated as B r . Hence the field H r can be related to B r defined by e rjm (B r ) =0. As a 
consequence of the Fermi statistics for the rapidities, £ r , m (A) < corresponds to occupied 
states (particles) and £ r , m (A) > to empty states (holes). Note that B r = oo for all r in 
the absence of external fields, i.e., H r = 0. 

In the limit T —>■ with fields H r the thermodynamics of the system is determined 
by £r,m(A) and with some tedious calculation, Eq. (|) for n = m can be rewritten as the 
following set of iV — 1 non-linearly coupled integral equations 

JV-l 

N 



er,m(A) + E^K ,ml * + e £ - /T ) = r(N M r) H r - 2tt JF r N (X), 



(10) 



s=l 



where Ki N ' m ^ is the Fourier transform of 



K[NM(,,\ - , m M sinh(mOT(r,g)f)sinh([Ar-maa;(r,g)]f) 

1 J sinh(mf ) sinh(iVf ) ' r ' s - 1 J 



The free energy per site is given by 



N-l 



f[N,m]( T> H) = f [N ' m] (0, 0) - £ T / F S N {\) ln(l + e e - m/T 



8=1 



(12) 



where 



(0,0) =—5: 



2J 



N 



k=\ 



k-1 

N 



with \l/ being the digamma function. 

We consider a level splitting of the iV-fold multiplets into a ground multiplet 
of degeneracy j and an excited (N — j)-fold multiplet due to a small external field Hj. 
This means that Bj is finite, but all other B r ^j = oo. Hence we can rewrite the integral 
equation fllCf ) for r = s = j as 



s jjm (X) + K^ M * Tln(l + e^' T ) = ^ j) Hj - 2ttJF^(X). 



(13) 



To evaluate the low-temperature specific heat coefficient in field, it is sufficient to expand 
£j jTn (A) to order T 2 , i.e., e^ m ~ ef^ + T 2 Sj m . Using the Sommerfeld formula to expand 
ln(l + e £j - m ^ T ) at low T, we have then the follwing integral equations for Sj m and £p m 



efi(X) + 2j^dXK^ m \X-X')e^{X') y 



r[N,m] 



j(N-j) 



Hj - 2nJF 1 N (X) 



71- 

~6~ 



(o) 



[K(X + Bj)+K(X-Bj)). 



The free energy is expressed by 



f[N,m](T } H) 



/(0,0) 



rfAFf (A)e °i(A) 





dE (0) 


I 3 


dX 



-1 



Ff(B j ) + 2l dXF" \\)ef n {\) 



B 



\ 

) 

! 



(14) 



(15) 
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The integral equations for e!j m were solved by Schlottmann |J for very small field yielding 
the susceptibility, i.e., Eq. (P and the following useful relations 



e — 



j(N-j)H jG i N ^\z%), 



47rJGf' mi) (0)sinf 



de (0) 

,j,m 



dX 



B, 



Q 



2nj(N-j)H 1 
(0) 



N 2Q{N,mj) 
1 



Q. 



(16) 

(17) 



m 



m 



2 \2irBj 



2 V 2 ^ , 



In 



m 



2nB„ 



+ ■ ■ 



Note that the parameter Bj tends to infinity as Hj — > as expected. 

To obtain the in-field specific heat, it is necessary to solve the integral equations 
for ef^ m in (TTJ|). Defining f(X) = e^(A + Bj), the integral equation becomes 



roc 

V(A) + / 


dX'[K(X - A') + #(A + A' + Bj)]if(X' 


7T 2 


de (0) 


-i 


6 


dA 


[jjf(A)+A'(A + 2S i )]. 





This equation can be solved iteratively for large Bj since the integral kernel K(B 



Hence substituting if ~ y>i + <f 2 + • • 
^x(A) + / dX'K(X- A'Vi(A') 



71" 



& (o) 

uc j,rn 



dX 



K(X), 



(f2(X)+ dX'K(X - A> 2 (A') = ^(-A-25,). 
Jo 

These equations are of the Wiener-Hopf type and can be solved by standard methods, 
yielding for A > 0, 



00 gL> 
— < 
-oo 2tt 





7T 2 


de {0) 






6 


cZA 





1 



G 



(N,m,j) 



UJ 



(19) 



where G 



(N,m,j) 



UJ. 



is defined in Eq. (|j). This solution does not have <£>2 contributions 
since <^2 is of order or higher and we are interested only in the leading contribution 
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to the in-field ^ N ' m K Collecting the results and inserting those into the free energy (p~5|) . 
we straightforwardly obtain that 



lim lim 7 

>0T-»0 • 



[N,m] 



12J [ + V N R n 
o 



For an arbitrary level splitting the linear superposition principle can be used as long as 
all Hj are small. In particular for a Zeeman splitting Hj = H, 



liml im-,W"l = ^-l) 



12J 



1 



N 



(20) 



Now let us consider the case when < T < J and Hj = 0. In this case we can 



calculate the entropy by the method used by Babujian and Filyov et. al. flSfl . As 
mentioned before, all rapidities obey Fermi statistics. This implies that the entropy per 
site of the model can be written as 

N-l oo oo 

5=EE d\[{a Stn + ln(a s ,„ + <rj n ) - a s , n lna s , n - tr* lnrrj ] , (21) 

8=1 n=l J - 00 

where the densities of occupied (u s , n ) and unoccupied ( a^ n ) states satisfy the integral 
equations (0). To calculate the entropy, we analyze the two sets of integral equations (|6|) 
and (|9]) which have the same integral kernels P^ s and driving terms as mentioned 
before. Substituting A — > A — A In ^ in Eqs. (H) and (|S|) and differenciating Eq.(|D|) 
with respect to A, we have the interesting relations between the energy potentials and 
the densities of states, i.e., 



CTr,n(A) 



N de r 



''-f(£r,n), 4(A) 



N ds r 



■[i - /Ml, 



(22) 



4vr 2 J rfA J v r '"" r ' nv y 4vr 2 J dA 
where /(x) = (1 + e X//T ) _1 . Note that these relations are valid only for sufficiently low T 
and [0, oo] in A-space is mapped on [— oo, oo] in A-space. For the convinience let us define 
E r>n (A) = s r ^ n (\ — A; In Tp-)/T and use the relation ( 22]) to rewrite the entropy ([H]) as 



S 



AIT N—l oo 

EE 

--1 n=l 



2ll 2 J 



dE s , n [f(TE Sin )\nf(TE s , n ) 
+(l-f(TE S!n ))ln{l-f(TE s>n ))\, 



(23) 



where E r>n satisfies the following integral equations 

N-l ^ H _,_ P B S ,„ + 1V1 -l_ P« S ,»-l' 



F fAl y pN , 1n (l + e*-"-»)(l + e j 

A r , n (A) - ^ ^ m (1 + eB ._ liB)(1 + eBi+liB , 



— Jsin(— n)e *> 5, 
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with the asymptotic condition lim^oo — ^ = 0. These equations imply that E r r 
E r>n (oo) and E™™ — E rtn (—oo). Hence the integral equations for maximum or mimimum 
solutions are reduced to the difference equations such that 



?max 
J r,n 



N-l 



El 



rpmin 



\- min{r,s) 1 + e »."+i)(l + e 

2_ ^ [A - max(r, s)\ In — gma ^ Wi — E mo 



8=1 

N-l 



El 



N 

y-, min(r, s 

8=1 

— oo. 



:i + e i 



0(1 + e J 



s + l,n 



A' 



■[iV — max(r, s)] In 



1 +e ^,n+l)(l 



-) 



1 + 6 * — 1 ' n 



1 + e ^+l,n) 



The solutions of these equations are of the form, 

(r + n)(iV — r + n) 



r,n<m 



rpmin 
r,n>m 



In 
In 
In 



r(iV — r) 



- 1 



sinf^vr) sin(^#%) 



-m+N ' 



m+N 



r f ir) sinf#^7r) 



- 1 



sin 



(24) 



rra+iV " / \rn+N 

(r + n — m) (AT — r + n — m) 



r(A^ — r) 

Note that the solutions for N = 2 with arbitrary m were obtained by Babujian [14 
who has treated the 577(2) invariant Heisenberg chain with arbitrary spin S. The above 
solutions ([24]) correctly reproduce those for N = 2 with m = 2S. 

Now we can express the entropy ([23]) as following, noticing that E, 

In x ln(i — x 



nmax 
J r,n 



Tpmin 
r,n+mi 



s 



AT r f 1 N—l m 

_iVi_ y y 

2-n 2 T ^ ^ Jn 

J r =l n=l J ° 



dx 



1 



+ 



(25) 



where 



sin 



m+N 



7r) sinf-^—^) 

' v m+N ' 



sinf^vr) sin(^=f±^7r) 



-m+N ' 



m+N 



We evaluate numerically the sums of the dilogarithmic functions, yielding 

7T 2 m(N 2 - 1) 



N-l m 

EE 

r=l n=l 



dx 



lnx ln(l — x) 
1 — X X 



3 m + N 



(26) 



Note that this sum rule is exact even though we have obtained it numerically. Hence the 
zero-field specific heat coefficient 



.. .. lNm] N(N 2 -l)m 
hm hm 7 lJV,mJ = V, 7 . 

oJ(m + A) 



which agrees with the result from the conformal field theory. It states that the linear 
specific heat coefficient 7 = HH, where the central charge c = m(N 2 — l)/(m + N) 



|T7| , |18| and the group velocity v = =z± |igj for SU(N) Heisenberg antiferromagnet. Note 
that for a sufficiently low temperature the all excitations have the same group velocities. 

The ^ N ' m ^ in finite field is less than the one in zero field as expected since the 
entropy is reduced in the ordered system. This result also reproduces the previous results 
for SU(2) invariant antiferromagnetic Heisenberg chain flID| , PH , noticing that the coupling 
constant J differs by a factor 2. 
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